Phase diagram of the disordered Bose-Hubbard model 
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We establish the phase diagram of the disordered three-dimensional Bose-Hubbard model at unity 
filling which has been controversial for many years. The theorem of inclusions, proven in Ref. 
states that the Bose glass phase always intervenes between the Mott insulating and superfluid phases. 
Here, we note that assumptions on which the theorem is based exclude phase transitions between 
gapped (Mott insulator) and gapless phases (Bose glass) . The apparent paradox is resolved through 
a unique mechanism: such transitions have to be of the Griffiths type when the vanishing of the 
gap at the critical point is due to a zero concentration of rare regions where extreme fluctuations of 
disorder mimic a regular gapless system. An exactly solvable random transverse field Ising model in 
one dimension is used to illustrate the point. A highly non-trivial overall shape of the phase diagram 
is revealed with the worm algorithm. The phase diagram features a long superfiuid finger at strong 
disorder and on-site interaction. Moreover, bosonic superfiuidity is extremely robust against disorder 
in a broad range of interaction parameters; it persists in random potentials nearly 50 (!) times larger 
than the particle half-bandwidth. Finally, we comment on the feasibility of obtaining this phase 
diagram in cold-atom experiments, which work with trapped systems at finite temperature. 

PACS numbers: 05.30.Jp, 63.50.-x, 03.75.Hh 
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I. INTRODUCTION 

The behavior of interacting bosons subject to static 
disorder is a fascinating subject whose study started more 
than 20 years ago 0, 0] ■ An important question raised in 
these papers is whether a direct transition between the 
gapped Mott insulating (MI) and superfluid (SF) phases 
is possible in the presence of disorder. Fisher et al. Q 
argued that a direct transition was unlikely, though not 
fundamentally impossible. Since then, the issue was a 
topic of hot debate with numerous analytical, computa- 
tional, and experime_ntal results reaching contradicting 
conclusions 

[i3,[ii,[il,[2c 

of direct 
mate approaches 
scenario! 

In Ref. [l| a final verdict was cast by proving ana- 
lytically that for any generic bounded disorder a direct 
transition between a superfiuid and a gapped insulating 
phase is not possible. Generic disorder is characterized 
by an arbitrary non- vanishing probability distribution of 
disordered fields within the bounds. Careful direct nu- 
merical simulations were in line with this prediction: in 
the presence of disorder, no matter how small, a Bose 
glass (BG) phase always intervenes between the super- 
fluid and Mott insulator phases. The Bose glass phase is 
an insulator with localized particle states at the chemical 
potential. Depending on system parameters these states 
can best be described either as localized single-particle 
levels or as isolated superfluid lakes. While the Bose 
glass docs not allow for phase coherence to extend over 



and some approxi 
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the entire system, it is characterized by a finite density of 
states and thus a finite compressibility and gapless par- 
ticle and hole excitations. The result of Ref. [l| comes 
as a simple corollary of the theorem of inclusions, which 
states that for any transition in a system with generic 
disorder one can always find rare regions of the compet- 
ing phase on either side of the transition line, provided 
the position of the line depends on the disorder distribu- 
tion function. However, there is a certain subtlety, if not 
a contradiction: The theorem seems to exclude any tran- 
sition between gapless and gapped phases in disordered 
systems, and the question arises of how to reconcile the 
theorem with the phase transition between the gapped 
Mott insulator and the gapless Bose glass phase. 

Previously it was conjectured 0, [H, [23|, but never 
proven rigorously, that the Mott insulator - Bose glass 
transition occurs when the bound A on disorder in the 
local chemical potential equals Eg/2. Here £"3/2 = 
min(£'p, i?/i) is the smaller of the particle [Ep) and hole 
{Eh) excitation gaps in the ideal Mott insulator (assum- 
ing that one works in the grand-canonical ensemble) ^ . If 
we denote by ^+ and ^_ the chemical potential thresh- 
olds for doping the Mott insulator with particles and 
holes respectively, then Ep = /i+ — /i, Eh = fJ- — It--- The 
gap for creating a particle- hole excitation (the MI gap), 
/i+ — is independent of the global 



E„ 



En + Eh 



^ In Ref. [23II . the convexity of free energy as a function of et — the 
crucial assumption in that paper — is a conjecture that might 
hold for the Bose-Hubbard model, but is incorrect in general, as 
shown by several counter examples. 



2 



chemical potential fi. At zero temperature, the chemical 
potential of the Mott insulator state with integer fill- 
ing factor can be anywhere between the two thresholds 
leading to an ambiguity in the value of Eg/2- The ambi- 
guity is absent in the canonical ensemble, where particle 
and hole excitations can be created only in pairs, to pre- 
serve the total number of particles. The grand-canonical 
counterpart of the canonical situation corresponds to the 
chemical potential being kept in the middle of the gap, 



H = — /x_)/2, in which case Ep = E^ 



E. 



gf2 



Eg/2. 



Therefore, below we always assume this choice of ii. 

The above-mentioned Ac = conjecture is based 

on the assumption that the state remains gapped for 
A < Eg/2- For A > Eg/2 the state can be shown to 
be gapless, because rare statistical fluctuations guaran- 
tee the existence of arbitrarily large homogeneous regions 
with disorder mimicking chemical potential shifts exceed- 
ing particle or hole gaps. In other words the conjecture 
was that the transition is of the Griffiths type. An al- 
ternative scenario would claim that the transition point 
happens at smaller values of A due to subtle interplay 
between disorder and interactions. 

In this paper, we show that the theorem of inclusions 
forces one to conclude that the Grifhths-type scenario 
is the only one possible for the gappcd-to-gaplcss tran- 
sitions. That is, the vanishing of the gap at the critical 
point is exclusively due to a zero concentration of rare re- 
gions in which extreme fluctuations of disorder reproduce 
a regular gapless system. In the vicinity of the critical 
point, the gapless phase must necessarily be "glassy" , be- 
cause it consists of large gapless (in our case superfluid) 
domains embedded in a gapped state. The absence of 
phase coherence between domains is caused by their di- 
verging distance between at the critical line. To illustrate 
these general conclusions, we consider the exactly solv- 
able random transverse field Ising model in one dimen- 
sion. 

Though the topology of the phase diagram for the 
Bose-Hubbard model is fixed by theorems, it is both in- 
teresting and important to determine transition lines and 
properties of phases numerically. In particular, this is 
necessary for revealing potential difficulties in observing 
and identifying the phases. To this end, we have cal- 
culated the full phase diagram of the disordered three- 
dimensional Bose-Hubbard model, shown in Fig. [Tl by 
quantum Monte Carlo simulations based on the worm 
algorithm [2^, This phase diagrams shows a few 

remarkable features: an infinite slope of the superfluid 
- Bose glass line Ac{U), in the weakly interacting gas 
U/t ^ 1, as predicted by the scenario of percolating su- 
perfluid lakes developed in Ref. [2^, and an enormous 
scale for the superfluid - Bose glass transition, A/t ^ 300 
at intermediate coupling strength, I ^ U/t ^ 30. Here 
U is the strength of the on-site repulsion between bosons 
and t is the amplitude of hopping transitions between 
the nearest neighbor sites (see Fig. [Ij. The percolation 
character of superfluidity in the vicinity of the superfluid 
to Bose glass transition, is most likely the reason for the 
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FIG. 1: Phase diagram of the disordered three dimensional 
Bose-Hubbard model at unity filling. In the absence of dis- 
order, the system undergoes a quantum phase transition be- 
tween SF and MI phases. The presence of disorder allows 
for a compressible, insulating BG phase, which always inter- 
venes between the MI and SF phases because of the theorem 
of inclusions [3]. The transition between MI and BG is of the 
Grifliths type, as an exception implied by the theorem. At 
— > 0, the SF-BG transition line has an infinite slope [25| . 



enormous scale. In this range of parameters, the localized 
states have a localization length of the order of one lat- 
tice spacing, as opposed to the picture of large superfluid 
lakes of Ref. [2|. 

The nature of the transitions and small superfluid frac- 
tion in the SF phase have profound implications for the 
experimental observation of the phase diagram. We focus 
here on cold-atom experiments, where recent experimen- 
tal claims are partly in line, partly in contradiction with 
the phase diagram shown above. We argue that present- 
day cold-atom experiments face numerous difficulties in 
obtaining the full phase diagram; for example, the Grif- 
flths type Bose glass ~ Mott insulator transition requires 
macroscopically large system sizes to properly identify 
the Bose glass phase. We also provide arguments why 
experiments seem to have missed the superfluid 'finger' 
above the Mott insulator in Fig.[Tl though the right scale 
for the transition between the superfluid phase and the 
Bose glass phase for very strong disorder has been re- 
vealed 

The paper is organized as follows. In Sec. [Til wc in- 
troduce the model and recapitulate the theorem of in- 
clusions. The transition between the Mott insulator and 
Bose glass phases is discussed in Sec. IIIII and illustrated 
by the exactly solvable random transverse Ising model 
in one dimension. We proceed with a discussion of the 
full phase diagram in Sec. IIVI and results of cold-atom 
experiments in Sec. [Vj The conclusions are presented in 
Scc.NB 
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II. MODEL AND THEOREM OF INCLUSIONS 

The disordered Bose- Hubbard model on a simple cubic 
lattice is defined the Hamiltonian 

H = -i^a]afc+^(ej - ^i) hj + ^^hj{nj-l), (1) 

(ok) 3 j 

where Oj is the creation operator of a boson on a site j; 
the symbol (. . .) denotes summation over nearest neigh- 
bor pairs of sites; fij = a^aj is the boson density opera- 
tor; and ej is the disordered on-site potential. Without 
loss of generality, we take ej to be independent random 
variables distributed according to the probability density 
p{e/A). The probability distribution satisfies the normal- 
ization condition J -^^ dup{u) = 1. has zero first moment 

^^-^ duup{u) = (otherwise it is absorbed in the defini- 
tion of and is taken to be bounded, that is p{u) — 
if |w| > 1. Formally, the disorder bound A and the disor- 
der distribution dispersion 5 are independent parameters. 
For the most common choice of the uniform distribution 
p{u) = const (used in our numerical simulations as well), 
we have (5^ = A complete characterization of dis- 

order is based on infinite number of parameters fixing the 
shape of p{u) . One may also add parameters which con- 
trol correlations between potentials on different lattice 
sites, etc. Collectively, we denote all these parameters 
by ^ and identify them with the definition of a particular 
model of disorder. 

Suppose now that a disordered system, described by 
the Hamiltonian similar to Eq. ([T]), undergoes a transi- 
tion from phase A to phase B — let us for the moment 
not specify the nature of phases — as the disorder bound 
A increases, and that the transition happens at a criti- 
cal point Ac, see Fig. [21 One obviously expects that Ac 
depends on the disorder strength 5, correlations between 
the sites, etc. For example, if correlations are long ranged 
and p{u) is very close to a (5-distribution, we hardly have 
any disorder at all no matter what the bound is, while for 
the uniform uncorrelated distribution disorder can radi- 
cally change system properties for large A. Let us define 
the notion of the generic A-B transition as a transition 
with some Ac(^) dependence. Figure [5] then proves that 
close to the transition, if A > Ac, there exist domains in 
phase B which locally look like phase A. Indeed, in the 
system described by Eq. ([T|) at A > Ac one can always 
find statistically rare domains where disorder realization 
is such that, within that domain, represent a typical 
realization of another disorder distribution with a bound 
< Ac, sec the dashed line going from B to A in Fig. [51 
Locally these domains are in the phase A. The proba- 
bility of observing such domains decreases exponentially 
with their size. 

What is more important is that close enough to the 
transition when A < Ac, there exist domains in phase 
A which locally are in phase B. At first glance, this is 
hard to justify, because the argument of the previous 




FIG. 2: A sketcli of tlie generic pliase transition line between 
some pliases A and B in the plane of disorder distribution 
parameters A and 5, wliere A is the bound and ^ is one of 
the infinite number of parameters characterizing the disorder 
distribution function, e.g. dispersion, and its spatial correla- 
tions. Dashed lines with arrows originate from points which 
determine disorder properties in the macroscopic (thermody- 
namic limit) system and end at points which characterize dis- 
order parameters in an arbitrarily large, but finite, domain as 
a result of a rare statistical fluctuation in the same system. 

paragraph can not be used as A has to be less than Ac 
everywhere. However, if we think in terms of all possible 
models of generic disorder we recognize that the actual 
value of Ac depends on the details of the distribution 
function p{u). This implies that it is always possible to 
choose ^ such that Ac(^) < A, and thus there are going to 
be domains in phase A, albeit exponentially rare as they 
get larger, where local disorder is indistinguishable from a 
typical realization of disorder with the distribution p^{u) 
and Ac(^) < A, see the dashed line going from phase A 
to phase B in Fig. [21 These domains will contain phase 
B. 

The above argument shows that it is not possible for 
phase A to be gapped if phase B is gapless. Indeed, the B- 
domains of arbitrarily large size within phase A guarantee 
that phase A is also gapless. As a consequence, no direct 
transition between the gapped Mott insulator phase and 
the gapless superfiuid phase is possible. 



III. GRIFFITHS TRANSITIONS 

A. An exception implied by the rule 

The theorem of inclusions rests on the dependence of 
the critical point on disorder properties such as its dis- 
persion, correlations, etc. Still one expects the gapful- 
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to-gaplcss MI-BG transition to exist, in apparent con- 
tradiction with the theorem! The paradox is resolved by 
considering the only remaining possibility, namely, that 
the transition point Ac does not depend on ^! In this 
case one cannot use arguments of the previous section 
to prove that in the vicinity of the transition point one 
can find arbitrary large domains of gaplcss phase B (we 
identify B with the Bose glass) inside the A (identified 
with the Mott insulator) ^. 

The transition which depends only on the bound A 
cannot be linked to any local physics, because as the 
dispersion 6 goes to zero the system becomes indistin- 
guishable from a pure one on larger and larger scales. 
This forces one to conclude that the transition mecha- 
nism itself is necessarily based on rare statistical fluctua- 
tions which explore the possibility of reaching the disor- 
der bound at all sites on larger and larger scales. Suppose 
that a gapped phase can be rendered gapless by applying 
a regular external field H. For the Mott insulator such 
a field is a global chemical potential shift Sfi; whenever 
fi-\~Sfi is above /i+ or below /i_ the system is doped with 
particles or holes and enters the superfluid state. The 
pathological insensitivity of the critical value Ac on ^ is 
natural for this scenario of rare regions in which the dis- 
order fluctuation is reproducing a regular pure system in 
an external field. When the disorder bound allows one 
to reach the critical value of the field, a transition oc- 
curs. We recognize that this mechanism is nothing but 
the conjectured Griffiths type MI-BG transition when 
the vanishing of the gap at the critical point is due to 
an infinitesimal concentration of rare regions in which 
the fluctuation of disorder mimics a homogeneous chem- 
ical potential shift In the general case it can be any 
regular external field whose amplitude scales with A. 

We thus conclude that that gapless-to-gapful transi- 
tions in disordered systems are possible if, and only if, 
they are of the Griffiths type and the transition line is 
fully determined by the properties of a pure system. In 
this case the disorder bound protects the gapped phase A 
(the Mott insulator) from having rare regions of phase B 
embedded in it. At the same time, when A is only slightly 
larger than Ac, then phase B appears to be identical to 
phase A locally except that it has rare, well-separated 
regions containing a gaplcss pure system. This means B 
cannot be superfluid, i.e. it is a glassy state. 



B. Illustration 

To illustrate the arguments presented above, we con- 
sider a disordered one-dimensional quantum model that 
shares some features with Eq. ([1]), but is exactly solv- 



^ Given an infinite number of continuous parameters determining 
disorder properties the probabiUty that any particular model of 
disorder go is a minimum is zero. 



able. It is closely related to the the two-dimensional 
classical Ising model with bonds whose strength depends 
randomly on their position in one spacial direction while 
being independent of the position in the second spatial 
direction. This model was flrst solved in Ref. [23] ■ It 
is characterized by a high temperature gapped param- 
agnetic phase, low temperature gapped ferromagnetic 
phase, and the intermediate Griffiths phase 30]. The 
Griffiths phase is akin to the Bose glass phase in model 
Eq. ID), while the gapped phases are similar to the Mott 
insulator. The nature of the transition between these 
phases can thus be clarified with the help of the exact 
solution (see also the discussion in Ref. [3l|). 

The model we consider here is a continuum limit in 
the second spatial direction, which we interpret as time. 
Then it is equivalent to the random transverse field one- 
dimensional Ising model, which wc can write as 

H = J2(ta^a^+,+h,a]). (2) 

j 

Here tij and cr| are Pauli matrices, acting on a j-th site 
of a linear chain, while hj are random independent vari- 
ables. The probability distribution p{h) is taken to be 
uniform on the h € [Hq — A, Hq + A] interval, where Hq 
and A are positive parameters such that ho > A. In 
principle, one could also add disorder to the Ising cou- 
pling t, however, this is not needed for the purpose of our 
illustration here. 

Model ^ is solved exactly by the Jordan- Wigner 
transformation [s^ . which became standard for these 
types of problems. Let us briefly review this method. 
With the notations cr+ = {a"" +iay)/2, a- = {a'=-iay)/2 
we introduce the Jordan- Wigner fermions 

flj = aj e ^kKi k fc ^ Q,j = e ^fe<i . (3) 

They satisfy the usual fermionic anticommutation rela- 
tions 

|a],afc|^ = {aj,afc}^ = 0. (4) 

In terms of these, the transverse field Ising model be- 
comes 

= 5Z - +t ^ [a] - cij^ [a]^^ + cij+i^ . 

j j 

(5) 

This Hamiltonian has the standard Bogoliubov form fa- 
miliar from the theory of superconductivity and can be 
rewritten as 

H = Y^{a] a,)n,Jl^), (6) 
jk ^ '^^ 

where 
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and I? is a matrix defined as 



(8) 



The problem now reduces to diagonalizing the real sym- 
metric matrix Ti.. To do so, it is convenient to perform 
first a unitary (actually, in this case, orthogonal) trans- 
formation defined as 



This gives 



V 
2?^ 



(9) 



(10) 



We recognize in Ti. a random one-dimcnsional Hamilto- 
nian in the BDI symmetry class, according to the classifi- 
cation scheme of Ref. |33j . This means that Ti. is real and 
that there exists a matrix E^^, in case of Eq. (jlOp given 
by 







such that 



-n. 



(11) 



(12) 



The arguments of Rcf. [33[ relate the peculiar proper- 
ties of the spectrum of the Hamiltonian (|10p which are 
discussed below to the existence of the symmetry ([T^. 

The problem defined by the Hamiltonian Eq. ^TU\i to- 
gether with Eq. ^ was solved exactly over thirty years 
ago in Ref. [3J| by the transfer matrix techniques, now 
standard in one dimensional disordered systems. In prin- 
ciple we could use this solution to extract all the infor- 
mation we need about the random transverse field Ising 
model Eq. ([2]) and its phase transitions. Yet the solution 
presented in Ref. [s^l is still relatively involved. To illus- 
trate the main features of the phase transitions, we can 
go to the continuum limit of Eq. (jlOp . In the continuum, 
the corresponding problem was solved in Rcf. • Their 
method is very simple and versatile, so we would like to 
put it to use here. 

The continuum limit in the Hamiltonian Ti. occurs close 
to the center of the band or to the momentum tt, when hj 
do not deviate much from some average value h^. More 
formally, we need to further transform the Hamiltonian 
by the unitary transformation given by 



U = 











(13) 



which keeps the structure of Ti intact but with the matrix 
T> now given by 



T>jk — hjSjk — t5j,k-i 



-t i5j,k-i - Sjk) + (hj - t) Sjk- 

(14) 



Now it is clear that the continuum limit of our problem 
is given by the same matrix Ti of Eq. (fTU|) with 



V = ta 



d 

dx 



V{x) 



(15) 



Here the continuum variable x is taken to be equal to j 
and a is the lattice spacing, while 



Vix) 



ta 



(16) 



In the continuum, V{x) can be thought of as a spatially 
random potential. 

Now the methods of Ref. [3^ (as adopted for this prob- 
lem in Ref. [1^) can be brought to bear on this problem. 
One result is that the spectrum of Ti is fully gapped if 
V{x) is everywhere positive or everywhere negative. Re- 
calling the definition of V{x) via hj and properties of the 
probability distribution p{h), this implies ho > t + A or 
ho < t — A. In other words, hj are either all greater than 
t or all smaller than t. 

Consider, for example, the case oiV>0. Suppose one 
increases A until regions appear where V{x) < (or hj < 
t). As soon as they appear, the spectrum of Ti becomes 
gapless. The density of states for positive energies E can 
be computed using the following construction. Take all 
regions where V{x) < 0. Consider the probability P{E) 
that 



\nE> dxV{x), 



(17) 



where the integration goes over one of the continuous 
intervals [xi,a::2] where V{x) < 0. Then the density of 
states is given by 



PiE) 



dP{E) 
dE 



(18) 



with a > — 1. We expect this probability to be exponen- 
tially small in ln£', or 



F(i;)-exp[(a + l)ln£;], 
where a is some number. Then 

p{E) - E". 



(19) 



(20) 



The derivation of Eqs. ^ and is given in 

Ref. [1^, and this completes the exact solution. 

We are now in a position to fully describe the transi- 
tion from a gapped paramagnetic phase with all hj > t 
to the gapless Griffiths phase as the disorder strength A 
is increased. For A > /iq — t, large rare regions appear 
where V[x) ^ hj — t is negative. We can use Eq. pT]) 
to calculate a, whose precise value depends on the prob- 
ability distribution p{h) but which in general is equal 
to some large number decreasing as A is increased past 
ho — t. Thus the low energy states which appear as A is 
increased above the threshold ho ^ t will be suppressed 
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by a power law, according to Eq. ((20|) . The Griffiths 
phase we obtained in this way is characterized by gapless 
excitations whose density is suppressed at low energy. 
Sometimes such a phase is referred to as a phase with a 
pseudogap (similar to a Mott glass phase arising in sys- 
tems with exact particle-hole symmetry and off-diagonal 
disorder [10, 113). 

We observe that the transition from the gapped param- 
agnetic phase to the gapless (but glassy) Griffiths phase 
proceeds exactly via the route described in this paper. 
When A < ho — t, no disorder, no matter what the de- 
tails of its distribution are, can create gapless states. The 
transition to the Griffiths phase occurs when disorder is 
just strong enough to create regions where gapless exci- 
tations can reside, because in this region an effective field 
h = ho — A can be made arbitrarily close to the critical 
value he = t. Wc note that the difference between the 
Bose-Hubbard model Eq. (JlJ and the random transverse 
field Ising model Eq. ^ lies in the fact that Eq. ^ , even 
in the fully clean (no disorder) regime, does not have a 
truly gapless phase, such as the superfluid in the Bose- 
Hubbard model. Yet the fact that Eq. ^ has a critical 
point in the absence of disorder is sufficient to create a 
glassy Griffiths phase with gapless excitations described 
by the power-law density of states. 

Identical arguments describe the transition from the 
gapped ferromagnetic phase to the Griffiths phase if ho < 
t as disorder strength A is increased past t ~ ho- 



IV. GLOBAL PHASE DIAGRAM 

In view of ongoing experimental activity to study the 
physics of interacting disordered bosons in optical lattice, 
and to connect the limit of strong interactions where dis- 
order competes with the physics of Mott insulators with 
the physics of localization in weakly interacting systems, 
we performed first-principles quantum Monte Carlo sim- 
ulations of the model ([1]). The results at unity filling 
are presented in Fig. [TJ As predicted by the theorem of 
inclusions, the superfluid and Mott insulator phases are 
always separated by the Bose glass phase at any A > 0. 

An interesting feature of the phase diagram from Fig. [1] 
is the reentrant nature of the Bose glass - superfluid 
transition if the interaction strength is increased at fixed 
disorder (and as long as disorder is not too strong to 
suppress the SF phase completely), confirming previous 
studies d, [111 . Fig. [3] shows this reentrant behavior on 
the more familiar (/i -I- 6t)/U vs t/U phase diagram. 
The dashed-dotted line represents a line of unit filling, 
(n) = 1, as the interaction strength is increased for fixed 
disorder strength. 

If disorder strength is increased, the Bose glass - super- 
fluid boundary line moves off to the right, and the entire 
line representing the Bose-Hubbard model at varying U 
may end up inside the Bose glass phase. 

Somewhat surprisingly, we find that the supcrfiuid 
phase extends to the region of very strong interactions 



([i+6t)/U 




FIG. 3: The phase diagram (/i + 6t)JU vs t/U of the Bose- 
Hubbard model, adapted from Ref. [3[. A dashed-dotted line 
represents the process whereby the interaction strength U is 
increased at a fixed disorder A in a system with unity filling 
factor, in the regime of the reentrant SF-BG transition. 

and disorder where both U and A are about two orders 
of magnitude larger then the hopping amplitude. This 
narrow finger-like region has fragile superfluid properties: 
being surrounded by the insulating state it has to have 
low transition temperatures to the normal state and small 
superfluid density at T = 0. Indeed, both quantities go 
to zero at the phase boundary. The superfluid transition 
temperature at A/t = 65 and U/t = 60 in the middle of 
the finger base is as low as Tc/t sa 0.37(5). Correspond- 
ingly, weak coherence properties (small condensate frac- 
tion) arc expected in the finger region. Moreover, they 
can be observed only on sufficiently large scale, because 
the correlation length diverges along the boundary. 

In another region of the phase diagram in Fig. [1] for 
U 0, the data clearly indicates an infinite derivative 
of the Ac{U) curve, in line with the prediction 

Ac cx U^/^ (21) 

of Ref. |25| . The results of Ref. [1^ were based on Gaus- 
sian random disorder, as opposed to a bound disorder dis- 
tribution discussed in this paper. However, we find that 
nevertheless their arguments remain qualitatively correct 
in our case too. Indeed, the nature of the Bose glass - 
superfiuid transition at weak interaction strength when 
U <^ t, A t, as discussed in Ref. [2^, is based on per- 
colation between localized states with energies ii^ << A. 
In this energy range, for states with large localization 
length, the Gaussian character of disorder fiuctuations is 
guaranteed by the central limit theorem. 

Quantitatively, we find that extremely large disorder is 
necessary to localize bosons even when interactions are 
relatively weak {U/t = 1 in terms of the lattice model 
parameters corresponds to the gas parameter naf ~ 10"'' 
for the continuous weakly interacting gas with the s-wave 
scattering length a^). 

At moderate interaction strength and A 3> i we ex- 
pect that the transition between the superfluid and the 
Bose glass phase is still driven by percolation. Because 
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of the imposed commensurability the mechanism differs 
from the conventional Anderson locaHzation argument 
which would predict a critical disorder of the order of the 
bandwidth. For strong disorder all single-particle states 
are localized with the localization length close to unity 
(in terms of the lattice constant) [3§]. The local (site) 
Hamiltonian 



Hioc = (e - fJ')n + ^n{n - 1). 



(22) 



can be used to determine the site occupation number as 



U/2 + n-e 
U ' 



(23) 



which is valid if n > or e < U/2 + ji. Otherwise, n = 0. 
As long as A ^ J7 the density n can be considered as a 
continuous function of ^ and there is no need to take into 
account that n can only be integer. The average density 
is now 



1 

2A 



A 



U 



-de. 



Setting (n) equal to unity leads to 



/i 



-U/2 



(24) 



(25) 



A site will be occupied if its disorder lies within the 
{-A,U/2 + n) = (-A,-A + 2\/c7A) interval. The corre- 
sponding probability is a/CT/A. If we assume that super- 
fluidity requires that occupied sites form a percolating 
cluster, then for a simple cubic 3D lattice with the per- 
colation threshold pc ~ 0.31 [iO] we find the transition 
line at 



E>1 
A ~ 10' 



(26) 



This estimate is in good quantitative agreement with 
the Monte Carlo results shown in Fig. [T] for intermedi- 
ate coupling before the Mott physics becomes important 
at U/t > Uc/t = 29.34(2). 

In turn, the assumption made above relies on the fact 
that moving a boson from one occupied site to another 
requires energy of the order of which is independent 
of disorder, while moving it to an empty site requires a 
much larger energy of the order of A. We note in pass- 
ing that while the percolation scenario drives the system 
towards the Bose glass - superfluid transition and thus 
defines, with a certain accuracy, the position of the crit- 
ical line, the criticality of the transition is most likely to 
be universal everywhere on the phase diagram. 



V. IMPLICATIONS FOR COLD-ATOM 
EXPERIMENTS 

Recently, experiments with ultracold gases have ad- 
dressed the disordered three-dimensional Bose-Hubbard 



Q. 




100 150 200 250 300 350 

A/t 



FIG. 4: Superfluid density as a function of disorder strength 
A at fixed interaction strength U/t — 30 for a system size 
L = 8 X 8 X 8 and inverse temperature pt = 10. The low value 
of the superfluid density shows the fragility of the superfluid 
finger in the large portion of the phase diagram Fig. [T] 



0) 

Q. 




100 120 140 



FIG. 5: Superfiuid density as a function of interaction 
strength U at fixed disorder strength A/t = 65 for a sys- 
tem size L = 8 X 8 X 8 and inverse temperature l3t — 10. 
The low value of the superfluid density is indicative of weak 
coherence properties and low transition temperature to the 
normal state in the 'finger' region. 



model [281 . l4l| . The random potential is generated using 
a fine-grained optical speckle field with correlation length 
comparable in size to the lattice spacing, but the disorder 
realization is usually kept fixed. The system is probed by 
looking at interference images giving access to the con- 
densate fraction, uq, provided the time-of-flight duration 
is sufficiently long and uq is large enough to be resolved in 
a trapped system. Transport properties are obtained by 
measuring the motion of the centre-of-mass of the atomic 
cloud immediately after an applied impulse [i^ . 

There are numerous considerations one has to keep in 
mind when trying to compare any experimental data to 
theoretical predictions for the homogeneous thermody- 
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namic system. The optical speckles not only introduce 
diagonal site-disorder, but also effect the on-site repul- 
sion strength and the hopping amplitude. In addition, 
there is a parabolic confinement trap rendering the sys- 
tem mesoscopic and inliomogeneous. This means that 
there is often a mixture of phases in the trap, such as 
the wedding cake structure where commensurate Mott 
domains are separated by liquid regions. Finally, exper- 
iments are done at low, but finite temperature. All of 
this complicates a direct comparison with the theory. It 
is however believed that the experiments can capture the 
phases and the transitions to some degree. 

The lack of a genuine compressibility measurement and 
a direct measurement of the gap make it difficult for cur- 
rent experiments to distinguish between the Mott insula- 
tor and Bose glass phases (they o nly distinguish between 
superfluid and insulating phases [2^). Moreover, the na- 
ture of the Griffiths transition prevents any experiment 
from direct observation of the transition line, because this 
would require astronomically large system sizes. As dis- 
cussed above, on short scales the Bose glass phase docs 
appear identical to the Mott insulator phase. Since dif- 
ferentiating between the Mott insulator and Bose glass 
phases in the neighborhood of the Uc/t point is not pos- 
sible experimentally, we will discuss here only the super- 
fluid - Bose glass transition. 

Experiments find that disorder can induce a superfluid- 
to-insulator transition, but they see no evidence for a 
disorder-generated insulator-to-superfluid transition, in 
apparent contradiction with Fig. [TJ At this point we 
recall, that superfluidity in the finger region is easily de- 
stroyed even by small finite temperature, because even 
at the base of the finger at A/t = 65 and U/t — 60 
the transition temperature is only Tc/t « 0.40(5) - such 
low temperatures were never reported in the literature 
for the Bose-Hubbard model in the strongly correlated 
regime. Furthermore coherence is weak even for T < Tc, 
see Fig. [5l It is likely that both effects are important in 
understanding why this region will be missed in the time 
of flight image. 

On the positive side, experiments do find the super- 
fluid - Bose glass transition for C//i = 25 and A/[/ = 10 
or A/t = 250. From our single site localization argument 
in combination with percolation it is clear that there is 
no fundamental problem in observing this transition ex- 
perimentally at sufficiently low temperature, because it is 
dominated by the short-range physics (the size of the su- 
perfluid region shrinks at finite temperature) . This find- 
ing gives the right order of magnitude answer when com- 
pared to the phase diagram shown in Fig.[TJ We note that 
a precise determination of the transition point requires 
sufficiently big samples that arc uniform in the middle 
(the Monte Carlo results had to be extrapolated to the 



thermodynamic limit since our answers drifted about ten 
percent for the lowest system sizes studied), and needs 
an accurate thermometry to study temperature effects. 
Our Monte Carlo results for the superfluid density as a 
function of disorder at fixed U/t = 30, see Fig. |4] indicate 
that Ug is severely depleted at large disorder, and thus 
transition temperatures in this region are small. 

Some of the difficulties discussed here arc specific to 
the Mott physics at integer filling factor. Away from 
commcnsurability, cold atom experiments can probably 
be successful in discerning the insulating glassy phase 
from the superfluid one. 



VI. CONCLUSIONS 

Summarizing, we have shown that Griffiths-type tran- 
sitions form a unique exception to the theorem of in- 
clusions. This immediately implies that all the gapful- 
to-gaplcss phase transitions in disordered systems are of 
the Griffiths type, and, correspondingly, close enough to 
the critical point, the structure of the gapless phase is 
what can generically be referred to as Griffiths glass: The 
system of distinct gapless domains containing a regular 
gapless system embedded into the gapped phase. This, 
in particular, proves the Griffiths nature of the Mott- 
insulator-to-Bose-glass phase transition. With the local 
shift of the chemical potential being the relevant field 
closing the MI gap, the critical line is given by the condi- 
tion 2 A = Eg^ where A is the bound of disorder and Eg 
is the particle-hole gap in the pure system. We have also 
considered a particular example of exactly solvable ran- 
dom transverse field Ising model which perfectly agrees 
with the established general picture. 

The full phase diagram has been presented in Fig. [T] 
and we have discussed the reason behind extraordinary 
stability of the superfluid phase against disorder and in- 
teractions. In combination with the analytical results in 
one dimension [l^ , numerical results in 2D [l^l , and the 
theorem of inclusions [l|, this study completes a com- 
prehensive description of the disordered Bose-Hubbard 
model at zero temperature in all physically relevant di- 
mensions. 

This work was supported by the National Science 
Foundation under grants PIIY-0653183 and DMR- 
0449521 and by the Swiss National Science Foundation. 
The authors acknowledge hospitality of the Aspen Center 
for Physics where this work has been initiated. We are 
grateful to I. Aleiner, J. Chalker, B. De Marco, E. Dem- 
ler and P. Weichman for stimulating discussions. Part of 
the simulations were performed on the Brutus cluster at 
ETH Zurich. 



[1] L. Pollet, N. Prokof'ev, B. V. Svistunov, and M. Troyer, [2] T. Giamarchi and H. Schulz, Phys. Rev. B 37, 325 (1988). 
arxiv:0903.3867. 



9 



M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 
J. K. Freericks and H. Monien, Phys. Rev. B 53, 2691 
(1996). 

R. T. Scalettar, G. G. Batrouni, and G. T. Zimanyi, 
Phys. Rev. Lett. 66, 3144 (1991). 

W. Krauth, N. Trivedi, and D. Ceperley, Phys. Rev. Lett. 
67, 2307 (1991). 

L. Zhang and M. Ma, Phys. Rev. B 45, 4855 (1992). 
K. G. Singh and D. S. Rokhsar, Phys. Rev. B 46, 3002 
(1992). 

M. Makivic, N. Trivedi, and S. Ullah, Phys. Rev. Lett. 
71, 2307 (1993). 

M. Wallin, E. S. So/rensen, S. M. Girvin, and A. P. 
Young, Phys. Rev. B 49, 12115 (1994). 
F. Pazmandi, G. Zimanyi, and R. Scalettar, Phys. Rev. 
Lett. 75, 1356 (1995). 

R. V. Pai, R. Pandit, H. R. Krishnamurthy, and S. Ra- 
masesha, Phys. Rev. Lett. 76, 2937 (1996). 
B. V. Svistunov, Phys. Rev. B 54, 16131 (1996). 
F. Pazmandi and G. T. Zimanyi, Phys. Rev. B 57, 5044 
(1998). 

J. Kisker and H. Rieger, Phys. Rev. B 55, R11981 (1997). 
L F. Herbut, Phys. Rev. Lett. 79, 3502 (1997). 
N. Trivedi, in Condensed Matter Theories Vol. 12 (Nova 
Science publishers, 1997), pp. 141-157. 
P. Sen, N. Trivedi, and D. M. Ceperley, Phys. Rev. Lett. 
86, 4092 (2001). 

J.-W. Lee, M.-C. Cha, and D. Kim, Phys. Rev. Lett. 87, 
247006 (2001). 

N. Prokof'ev and B. Svistunov, Phys. Rev. Lett. 92, 
015703 (2004). 

J. Wu and P. Phillips, Physical Re- 
view B (Condensed Matter and Materials 
Physics) 78, 014515 (pages 11) (2008), U RL 
|http : //link. aps . org/ abstract /PRB/v78/eO 145 15 [ 
[22] U. Bissbort and W. Hofstetter, EPL (Euro- 
pliysics Letters) 86, 50007 (6pp) (2009 ), URL 
|http : //stacks . iop . org/0295-5075/86/50007| 



[23] P. B. Weichman and R. Mukhopadhyay, Phys. Rev. B 

77, 214516 (pages 39) (2008). 
[24] P. B. Weichman, Mod. Phys. Lett. B 22, 2623 (2008). 
[25] G. M. Falco, T. Nattermann, and V. L. Pokrovksy (2008), 

arXiv:0811.1269. 
[26] N. V. Prokof'ev, B. V. Svistunov, and I. S. Tupitsyn, 

Journal of Experimental and Theoretical Physics 87, 310 

(1998). 

[27] L. PoUet, K. Van Houcke, and S. M. A. Rombouts, Jour- 
nal of Computational Physics 225, 2249 (2007). 

[28] M. Pasienski, D. McKay, M. White, and B. DeMarco, 
arXiv:0908.1182v3. 

[29] B. M. McCoy and T. T. Wu, Phys. Rev. Lett. 21, 549 
(1968). 

[30] R. B. Griffiths, Phys. Rev. Lett. 23, 17 (1969). 

[31] D. S. Fisher, Phys. Rev. B 51, 6411 (1995). 

[32] T. D. Schuhz, D. C. Mattis, and E. H. Lieb, Rev. Mod. 

Phys. 36, 856 (1964). 
[33] M. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 
[34] T. P. Eggarter and R. Riedinger, Phys. Rev. B 18, 569 

(1978). 

[35] A. Comtet, J. Desbois, and C. Monthus, Ann. Phys. 

(N.Y.) 239, 312 (1995). 
[36] V. Gurarie and J. T. Chalker, Phys. Rev. B 68, 134207 

(2003). 

[37] T. Giamarchi, P. Le Doussal, and E. Orignac, Phys. Rev. 

B 64, 245119 (2001). 
[38] W. Krauth and N. Trivedi, EPL (Eu- 

rophysics Lett ers) 14, 627 (1991 ), URL 

http: //stacks . iop . org/0295^5075/14/627 
[39] B. Bulka, M. Schrcibcr, and B. Kramer, Z. Phys. B Cond. 

Mat. 66, 21 (1987). 
[40] M. B. Isichenko, Rev. Mod. Phys. 64, 961 (1992). 
[41] M. White, M. Pasienski, D. McKay, S. Q. Zhou, 

D. Ceperley, and B. DeMarco, Phys. Rev. Lett. 102, 

055301 (2009). 

[42] D. McKay, M. White, M. Pasienski, and B. DeMarco, 
Nature 453, 76 (2008). 



